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Abstract 



We calculate the mean number of subjets in quark and gluon jets in the 
final state of e + e~ annihilation. Since "quark" and "gluon" jets are scheme- 
dependent objects, we stress the importance of using the same definition as 
in experimental analyses. We define the jets using the k±_ algorithm at a 
coarse scale y±, and the subjets using a finer scale y , in the same algorithm. 
Gluon jets are anti-tagged by the presence of heavy quarks in both other 
jets. Our result is exact to leading order in a g , and resums leading and 
next-to-leading logarithmic terms in the ratio yi/yo to all orders in a s . 
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1 Introduction 



Leading logarithmic QCD predicts that the multiplicity of hadrons in a gluon jet should 
be larger than that in a quark jet by the ratio of their colour charges, C A /C F = 9/4. Thus 
one would naively expect the multiplicity in three-jet e + e~ events to be (2C F +C A )/2C F = 
17/8 times higher than in two-jet events. Experimentally, however this ratio is close to 
unity 0. The sub- leading corrections of relative size y/a s and even a s f2|, are known for 
the case of jets produced by a colour singlet source. Although they reduce this ratio, they 
are small at the energies of modern experiments, and certainly not capable of bringing 
it down to unity. When the jets are produced in a three-jet ensemble, rather than by a 
colour singlet, the equivalent corrections can be much larger 0], but without assumptions 
about the non-perturbative hadronization process, it is impossible to calculate their size. 

Considerable light was shed on this problem by the paper of Catani et. ai.Q, who 
considered an analogous problem that can be solved purely perturbatively, namely the 
multiplicity of subjets in a jet. A jet algorithm is run twice with different cutoff scales, 
Q and Qi, and the number of jets found at each, M. Q1 . A4o is then called the num- 
ber of subjets in an A^i-jet event. By considering the region A QCD <C Qo <C Q±, they 
were able to study the region in which logarithmic terms dominate, as in the hadron 
multiplicity problem, but which is fully under perturbative control. They found that 
colour coherence during the formation of the three-jet state was responsible for a deple- 
tion in the multiplicity in three-jet events relative to the naive, incoherent, expectation. 
Experimental results || have shown good agreement with the predictions of @] for not 
too small Qo, although there is clear evidence of an additional suppression in three-jet 
events as one approaches the non-perturbative region. 

More recently, using their excellent ability to tag heavy quark jets, experiments at 
LEP have been able to measure the subjet multiplicity in individual jets that are anti- 
tagged as the gluon jet in three-jet events||. This allows a more direct measurement 
of the multiplicity ratio between gluon and quark jets in three-jet events. Although 
the results were in good agreement with Monte Carlo event generators that incorporate 
coherence J7L |8], [|, it was not possible to compare with the results of M for two reasons. 
Firstly, their method is not able to keep track of which jet the registered subjets are 
part of, it simply counts the total number of subjets. Secondly, the experiments found 
that to obtain a sufficient lever-arm in Qi/Qo for the logarithmic behaviour to become 
important before the non-perturbative behaviour, Qi had to be rather large ~ Q, where 
Q is the total e + e~ annihilation energy, while M work in the Qi <C Q limit. 

In this paper, we correct both deficiencies by numerically integrating an expression 
that exactly reproduces both the leading order expression in a s and the next-to-leading 
order resummation of logarithms of Qi/Qo, but not logarithms of Q/Q\. Our result is 
therefore uniformly reliable for all Qi/Qo, and not too small Q\. In section 2, we give 
details of the event definition, which is matched as closely as possible to the experimen- 
tal one. In section 3, we discuss the logarithmic behaviour and define the numerical 
treatment, which allows us to resum logarithms of Qi/Qo, without losing the exact 
treatment of the Q\ ~ Q region. We compare our results with the data of |5| and find 
good agreement where the hadronization corrections are expected to be small. Finally in 
section 4 we briefly discuss the dependence on the energies of the individual jets within 
the three-jet sample. 



2 The Event Definition 



In order to be able resum logarithms of the jet resolution scale, jets must be defined 



using an algorithm such as the k± algorithm fllOf in which the phase-space for multiple 
emission factorizes in the same way as the matrix elements. This allows jet multiplicities 
to be resummed to at least next-to- leading logarithmic accuracy to all orders in a s |4], [Ljfl . 

The algorithm defines for every pair of particles in the final state a 'closeness' mea- 
sure, 

_ fcjjj _ 2min(E i ,E J -) 2 (l -cosgfj) 

Vij- q 2 - Q2 ■ UJ 

If the smallest value of y^ is below a cutoff y cut , the pair i,j is merged into a single 
pseudoparticle with momentum 

p ij =Pi+p j , (2) 

and the process is repeated, with pseudoparticles treated on an equal footing with 
particles. If all y^ are above y cu t the algorithm is terminated and all remaining particles 
and pseudoparticles are called jets. 

To define subjets, the algorithm is first run with a cutoff y cut = yo, and the resulting 
particles and pseudoparticles are called subjets. The algorithm is then continued with 
a cutoff y cnt = yi > yo, and the resulting particles and pseudoparticles are called jets. 
During this second stage of clustering, the fate of each subjet is kept track of, giving 
the number of subjets 'inside' each jet. 

At present, the exact matrix elements for the final state of e + e~ annihilation are only 
known up to 0(a;?)|112|. Since the multiplicity of subjets in an n-jet event only becomes 
non- trivial at 0(a^), we are confined to n < 3 in exact calculations. Furthermore, 
since the two-jet final state does not contain gluon jets, we are only interested in n > 3. 
Therefore, we can only calculate a single non-trivial term, the leading order expression 
for the subjet multiplicity in three-jet events. This can be easily found by numerically 
integrating the exact matrix element within the phase-space region given by the event 
definition. 

So far, we have fully defined the subjet structure of the event, but to obtain the 
maximum amount of information from such events, we must also define the flavour 
structure. Namely which jets are quark jets and which are gluon jets. Beyond leading 
logarithmic order, the results will depend on this definition, so we must use the same 
definition as in experimental analyses. This is done using heavy flavour tagging: in 
three-jet events in which two of the jets contain heavy quarks, the third is called a gluon 
jet. Then, using the subjet multiplicity in three-jet events, Ai 3 , and in gluon jets, M. g , 
the subjet multiplicity in quark jets can be found by assuming that every three-jet event 
consists of two quark jets and a gluon jet, 

M q = \{M*-M B ). (3) 

Using this definition, M. q is a physical quantity, since both and M. g are. 

At there are two sub-processes to consider: qqgg and qqq (/) q (/) . For each case, 

one must consider what flavour to call a jet consisting of a pair of partons of given 
flavours. When one of the partons is a gluon, this is straightforward: the flavour of the 
jet is the flavour of the other parton. However if both partons are quarks or antiquarks, 
one must carefully consider how the flavour-tagging would assign the jet, recalling that 



a gluon jet is defined as the untagged one in an event in which both the other jets are 
tagged as quark jets. 

In the case of qqgg, the combination q + q should be called a quark jet, since neither 
of the others is a quark. 

In the case of the qqq (/) q ( ' ) final state, q + q is called a gluon jet, since both the others 
are quarks. For the combinations q + q (/) and q + q', the outcome will depend on the 
exact experimental procedure and the corrections that are applied to it. We make the 
simplifying assumption that each jet is equally likely to be tagged, so these are shared 
between gluon and quark jets in the ratio 1 to 2. These contributions make such a small 
contribution to the total cross-section that varying this assumption does not affect the 
final result significantly. 

It is straightforward to integrate the leading order matrix element according to these 
definitions, and we obtain the results shown later, in figure [1|. 



3 Resumming Large Logarithmic Terms 

When the ratio of cutoffs, yi/yo, becomes large, double-logarithmic terms appear in the 
perturbative expansion, M. ~ a™ log 2n Vi/vq- Such terms must resummed to all orders in 
a s to give a reliable prediction. Although this was done in [|J, they used the simplifying 
limit yi <C 1, and were not able to keep track of which jet contained which subjet. In 
this section we correct both deficiencies, by using a very similar trick to that first used 
in (T^j for the subjet multiplicity in hadron collisions. 

To illustrate the general method, we first work in logarithmic approximation, valid 
for We use the usual notation L 0) i = — log?/o,i- A three-jet configuration 

in which the gluon makes an angle 9 with the nearer quark can increase the subjet 
multiplicity in five ways: gluon emission from the opposite quark, from the nearer quark 
at angles larger than 9, angles smaller than 9, emission from the gluon, and gluon 
splitting to quarks. To next-to-leading logarithmic accuracy, one can neglect recoils in 
the emission, and the multiplicity can be easily found, 

M 3 -3 = la s (L -L 1 ){C F (L + L 1 -3)+C F (lL 1 -l) + C F (L + lL 1 -2) 

+C A (L -lL 1 - 1 4) + lN f }, (4) 

where a s = a s /2ir. The separate terms correspond to the different contributions just 
mentioned, and add up to the expression of j|. 

Since there is one set of terms proportional to the colour charge of a quark, and 
another to the colour charge of a gluon, it is tempting to assume that these are the subjet 
multiplicities in the quark and gluon jets respectively. However, simple kinematics shows 
that this is not the case. For example, a gluon emitted from the nearer quark at a larger 
angle than 9 could be merged with either the gluon or the quark, depending on the 
relative azimuths. There are non-trivial azimuthal correlations in the soft limit of the 
qqgg matrix element, which are not retained in the usual evolution equations because 
they average to zero. Therefore the resummed result of M cannot predict the number 
of subjets that are merged into each jet, although it can predict the total number of 
them. We shall shortly see how this problem can be solved without having to define 
new azimuth-dependent evolution equations. 



Similarly, emission from the gluon can be merged with the quarks, depending on 
the azimuth and opening angle. However, simple kinematics also shows that emission 
from any jet at angles smaller than 9/2 will always be merged with that jet. We make 
use of this fact by introducing an arbitrary parameter \i that separates the intrajet and 
interjet regions. We consider \i to be fixed and in the range Qo < n < Qi/2f], although in 
general one could make it a function of the three-jet kinematics. We then call emission 
from a jet of energy E at angles smaller than fi/E intrajet emission (note that 9 must 
be larger than Qi/E), and at angles larger than \i/E interjet emission. This gives us 
multiplicities 

■Mntra — 3 = ia s {2C F (L 2 — + C A (L^ — y£ M o) + §^/-^M)}> (5) 

where L^ = log fi 2 /Ql, and 

AT intcr = |a s {2C F [L M0 (2L 1 + 2L Ml ) + L /il (2L 1 + L„ x - 3)] 

+C A [L M0 (i£i + 2L Ml - 1) + V(§Li + V - f )] + fiV/V}, (6) 

where L^i = \ogQl/fi 2 . We make the distinction between results for physical multi- 
plicities, which we denote A4, and theoretical multiplicities, which we denote N '. Note 
that 

A^ntcr + AT intra = M, (7) 

with no dependence on //. Since \x is an entirely arbitrary parameter separating different 
treatments of the same phenomena, it should cancel in all physical quantities, at least 
to the accuracy to which we are working. However, in general it will not exactly cancel, 
so varying it will give a rough indication of the size of neglected terms, and hence 
the accuracy of our calculation. In this sense, // plays a similar role to the familiar 
renormalization and factorization scales. 

We turn now to the problem of how to resum the large logarithms of Qx/Qo to 
all orders in a s , while keeping exact control of the three-jet kinematics, ie. without 
assuming Qi <C Q. From equations (|5|) and (||), it is clear that the only dependence on 
Qo is through L^ , so we aim to resum all terms of the form a^L 2 ^ (leading) and a™L 2 /^~ l 
(next-to-leading)[|. We do this separately for the intrajet and interjet components. 

Resumming the intrajet multiplicity 

The structure seen in equation ([|) continues to all orders, with the intrajet multiplicity 
contributing two logarithms of /i 2 /Qo f° r ever y power of a s . Therefore in resumming 
it, we must keep track of next-to- leading corrections to the evolution. By explicitly 
setting up the evolution equations, one finds that the intrajet multiplicities are identical 
to those in jets of the same flavour formed from a colour singlet at scale fi resolved 
with a transverse momentum cutoff Qq. Results for these quantities, which we denote 
NqiQoiH) and Af 9 (Qo, /i) for quark and gluon jets, were given to next-to-leading loga- 
rithmic accuracy in jLlJ. We do not bother repeating the definitions here, but it is worth 

*If Qo > Qi/2, there are no large logarithms to resum, so we do not use the following procedure. In 
that case the leading order matrix element is perfectly sufficient on its own. 
TWe are implicitly assuming fi ~ Qi. 



noting their expansions in the threshold region (formally defined as the region where 
ot s L 2 Q is small but L^ is still large) 



Af q (Q ,fi) = l + as^fiCpLjo-fa^oJ+Ofe 2 ), (8) 
N b (Qq,h) = l + a s (fi)[±C A Ll -lbL^ ]+O(a 2 s ), (9) 

where b = ^C A — In terms of these multiplicities, the total intrajet multiplicity is 
simply given by 

A/" intr a = 2A/",(Q , fl) + N g {Q Q , fl). (10) 

Note that \i naturally arises as the scale to use for a s when we evaluate the fixed-order 
expression. 

Resumming the interjet multiplicity 

Physically, there are two ways that the interjet multiplicity can be increased beyond 
leading order: either by multiple emission from the three-jet event, or by emission from 
the gluon emitted at first order. Since the leading order expression is single-logarithmic 
in L M o, raising it to any power results in terms that are negligible to next-to- leading 
accuracy. This means that multiple emission is negligible, and to retain next-to-leading 
accuracy overall we only need to keep track of leading logarithms arising from emission 
from the leading order gluon. 

As we already mentioned, to keep track of which jet an interjet subjet was merged 
into, we need to include azimuthal correlations between different emissions. Thus it 
seems that azimuth-dependent evolution equations are needed. However this is not 
the case, because to leading logarithmic accuracy, the additional partons inside the 
interjet subjet can be considered collinear, and so are always combined together by 
the jet algorithm before being combined en masse with one of the main jets of the 
event. Therefore we only need to calculate the distribution of subjets between jets to 
leading order in a s , to give an expression to next-to-leading logarithmic accuracy to all 
orders in a s . This is identical to the trick we used in [|13|], where we calculated interjet 
terms (in that case generated by initial-state radiation) for the number of subjets in a 
hadron-collision jet, despite being unable to obtain analytical results for the kinematic 
configurations that generated them. 

In calculating the number of subjets in the interjet region, one encounters integrals 
of the general form 

A/inter ~ I dk\ 0^, $)JV fl (Q , (11) 

where k± is the transverse momentum of the interjet gluon, $ represents the rest of the 
phase-space, namely an energy fraction and an azimuth, and G represents the (arbitrarily 
complicated) phase-space constraints. This can always be manipulated into the form 

A/i n ter = Antcr / -TrOC 8 (k±)J\f g (Q , k±) + a s (n)B inter , (12) 



where Ai nter and -Bmter are arbitrarily complicated non-logarithmic functions, and we 
have neglected M g in the second term, which is justified to next-to- leading logarithmic 



accuracy. Noting that the leading order expansion of M g is just unity, we obtain the 
exact result to leading order in a B , 

■Winter = «s()") Anter^O + «s G") Anter • (13) 

To calculate the integral retaining J\f g , we use the evolution equation to leading loga- 
rithmic order, 

A^Qo, fj) = C A l -j±a s (k ± )Af g (Q , k ± ) (14) 
where the prime denotes logarithmic differentiation with respect to the second index| 

M g (Qo,Q) = Q 2 ^K(Qo,Q). (15) 

Using this to rewrite (0), we arrive at the same result exact to both leading order in 
a s , and next-to-leading order in L^ for all orders in a s , 

A/inter = ^4^— M g {Q , /i) + « s (/i) -Sinter- (16) 



The final point concerns the evaluation of A inter and B inter . As we have stressed, they 
are arbitrarily complicated in general, and not amenable to analytic evaluation. How- 
ever, as in flTSfl , it is possible to evaluate A/inter to next-to- leading logarithmic accuracy, 
without knowing their explicit forms. We note that the threshold expansion of Af g is 

N' g {Qo,n) = C A a s (fi)L, + O(a 2 s Ll ), (17) 
and that -Bi n t er is non-logarithmic, so the relation 

Os-D mter —— — — = a s i3 inte r [i-O) 



holds to the appropriate accuracy. Comparing equations ( |l3l) and (|16|), we arrive at the 
final answer for the interjet multiplicity in three-jet events, exact to leading order in a s 
and resumming next-to-leading logarithms in L^ to all orders in a s , 

■/ V inter — J V i n t er „ _ , \ , • {W) 

Put simply, this means that the full result is found by numerically integrating the 
exact leading order matrix element, and multiplying the result by the parton multipli- 
cation factor, 

jVgg^ , (20) 
C a q; s (//)L m o 

Explicit forms for A in ter and -Binter are never needed. It is worth stressing again, that 
this correctly keeps track of which jet each subjet is merged into to all orders in a s , in- 
cluding azimuthal correlations. Note that the parton multiplication factor has threshold 
expansion 

r a - Q ( \f = 1 + i C ^H + 0(a% (21) 

so rises from unity with increasing L^q. 

*Note that this is a factor of two smaller than the J\f' g we defined in |]l3| . 



The full multiplicity 



In combining the results for the two regions, we have the opportunity to make a final 
embellishment that avoids the need to differentiate between the interjet and intrajet re- 
gions when numerically evaluating the leading order matrix element. This is because the 
intrajet region gives a contribution of the logarithmic terms (|5]) plus a non-logarithmic 
term. Multiplying non-logarithic terms by the parton multiplication factor gives negli- 
gible next-to-next-to-leading logarithmic corrections. Therefore the full expression for 
the multiplicity reads 



t-2 3^f T 



M = a s {fi)\^M as -2 

+2Af q {Q ,fi)+Ar g {Q ,fi), 



A^(Qo^) 
C A a s (fx)L 



(22) 



where all the dependence on \i has been explicitly displayed. 

As a final check of this formula, we show that it is independent of fi to the appropriate 
accuracy by taking the logarithmic derivative, 



AM 
djj? 




CaL^q — 



A^(Qo,/i) 
C A a s (fi)L 



\C A L 2 ^ — \bL^ 



x 



<(Q ,/i) A/-'(Q ,/i) 



+2Af'(Q , f i)+Af'(Q ,ti). 



(23) 



The first and third terms cancel each other exactly to leading order in a s and, using the 
leading logarithmic result 

A/;(Qo,A0 = ^Wo,A*), ( 24 ) 

to next-to- leading logarithmic accuracy for all orders in a s . The two terms in the 
second bracket of the second term exactly cancel to leading order in a s , but at first 
sight it appears that there are non-zero terms at next-to-leading logarithmic accuracy 
(a s (/i) x L^qX a next-to- leading logarithmic function). However, note that Ai as contains 
Lq terms with the same coefficients as the functions, so for Qq fi ~ Qi the 
first bracket is actual single logarithmic, and the whole term is zero to next-to-leading 
logarithmic accuracy. 

Therefore the expression ( |2"2"D is indeed independent of fi, as we claimed, and gives the 
full result exact to leading order in a s , resumming leading and next-to- leading logarithms 
to all orders in a s , and correctly keeping track of which jet the subjets get merged into. 

In figure [l]a we show Eq. (^) in comparison with the data of for y\ = 0.1. We fix 
Aqcd = 250 MeV, and use the two-loop running coupling (although this is not justified 
at our accuracy) corresponding to a s (M z ) = 0.120. We see that agreement is good down 
to about y Q = 2 x 10~ 3 , which corresponds to transverse momenta of about 4 GeV. Below 
this value, the data rise above the predictions. The barely distinguishable pairs of curves 
are the result of varying the separation scale /i, with the natural choice, Qi/2, giving 



ft 




y„ 1=0.1) y (y,=0.1) 



Figure 1: TTie sub jet multiplicity in quark and gluon jets, with jet definition scale 
Hi = 0.1, in comparison with ALEPH data^. (a) Absolute values for gluon jets (solid) 
and quark jets (dashed) with \i = ^Jy{Q/2 and \i = ^/y{Q/A. Also shown are the leading 
order results (dotted), (b) The ratio of gluon to quark jets, with fi = ^Jy{Q/2 (solid), 
fi = y/yiQ/4 (dashed) and at leading order (dotted). 




the lower curve, and Qi/4 giving the upper. Also shown are the leading order curves, 
which fall below the data well before the full results. 

In figure |l|b we show the ratio of the gluon and quark multiplicities. It can be 
immediately seen that in the region in which the multiplicities rise above our predictions, 
yo < 2 x 10 -3 , the quark multiplicity is increased more than the gluon, suppressing the 
multiplicity ratio. The effect of varying the separation scale is more pronounced simply 
because the vertical axis is more compressed. As shown in |J, the disagreement between 
our results and the data can be attributed to hadronization corrections both at small 
yo and the intermediate yo values where the data lie slightly above our curve. In fact 
the parton-level predictions of the JETSET program shown there agree almost perfectly 
with our curves, while the hadron-level predictions agree almost perfectly with the data. 

In figure |2| we show the ^-dependence of the results, which is considerable, although 
the general features remain unchanged as yi is varied. If a measurement of the in- 
dependence could be made, it could be used to separate the interjet and intrajet com- 
ponents of multiplicity, without needing to explicitly measure the subjet directions. 
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Figure 3: As in figure [| but for jets in the energy range 30 < £j e t/GeV < 35. 

4 Jet Energy Dependence 

In the full results there is clearly a reduction in the difference between gluon and quark 
jets, relative to the leading logarithmic expectation. One might suspect that at least 
part of this is because the gluon jet tends to be the least energetic, despite the fact that 
the selection yi = 0.1 requires all three jets to have fairly similar energy. This can easily 
be checked by binning the events according to jet energy. 

Checking Eq. fl22f ) critically, one finds that it holds equally well for the multiplicity 
as a function of the three-jet kinematics. Specifically, none of the resummed components 
depend on the three-jet kinematics, so if A4 as is binned in jet energy, then M. will also 
be binned correctly. 

As an example, we show in figure |3] the results for jets with energy in the range 
30 < Sjet/GeV < 35. We see that the multiplicity ratio is indeed enhanced, particularly 
at larger y values. 

5 Summary 

By combining numerical integration of the leading order matrix element with analytical 
resummation of large logarithmic terms, we have obtained an expression for the subjet 
multipicity in 3-jet e + e~ events ( p2| ) that is uniformly reliable for all values of yo/yi, 
and correctly incorporates azimuthal correlations between jet and subjet directions. 
Where the hadronization corrections are estimated to be small the results are in good 
agreement with data. It is possible that modifications to the jet definition could reduce 
the size of these corrections, and the scale at which they begin to become important. 
The method used in this paper could easily be adapted to any definition in which the 
systematic resummation of large logarithmic terms to all orders is possible. It would 
also be straightforward to adapt it for any other three-jet selection or flavour-tagging 
method, since only the fixed-order matrix element depends on the three-jet kinematics 
and this is integrated numerically. 
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